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The ; < : I : : n i* 1 1 <>;if m« nr of interact ive. realistically varying aerosol in a global 
opcmafiona! forecasting system, as opposed to prescribed (fixed or * lini.it . : 
logically varying) aerosols, is a very difficult challenge that only recently Ire- 
gins to be addressed. Experiment results from a recent version of the NASA 
G EOS-5 forecasting system, inclusive of interactive aerosol treatment, are 
presented in this work. Four sets of 30 5-day forecasts are initialized from 
a high quality set of analyses previously produced and documented, to cover 
the period from 15 August to It) September 2006, which corresponds to the 
NASA African Monsoon Multidisciplinary Analysis (NAM.VJ A ) observing cam- 
paign. The four forecast: sets are at two different horizontal resolutions and 
with and without interactive aerosol treatment. The net impact- of aerosol, 
at times in which there is a strong dust outbreak, is a temperature increase 
at the dust level., and decrease in the near-surface levels, in complete agree- 
ment with previous observational and modeling studies. Moreover, forecasts 
in which interactive aerosols are included depict an African Easterly Jet at 
slightly higher elevation, and slightly displaced northward, with respect to 
the forecasts in which aerosols are not included. The shift in the AEJ po- 
sition goes in the direction of observations and agrees with previous results. 
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i. Introduction 

The role of the Saharan Air Layer (SAL) has been .suggested to be relevant to weather 
. forecasting over the tropical Atlantic since at least, t he early 70s [Carlson and Prospers). 

1072;. In addition to the thermal effect due to its intrinsically high heat content, the 
» role of dust, in terms of direct radiative effect lias been intensely investigated for several 
j', decades (e.g., Carlson and Benjamin, [1980]). Surface- and space- based measurements of 
» atmospheric: optical thickness (AOT) go back to the late 70s (e.g. [Carlson and Wendliug. 
» 1977]) but the possibilities of merging modern-era satellites with data such as the ones 

* coming from Moderate Resolution Imaging Spectrum! iometer (MODIS) or Cloud-Aerosol 
» Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) give a much more aceu- 

* rate understanding of aerosols' optical proper t ies (e.g. [Remer et. al. 2008; Omar et a!. 

« 2009]) 

» Among the various effects of dust over the tropical atmosphere, Dunion and Velden 
» [2003] suggested a role of dust unfavorable to tropical cyclogenesis. One point of their 
.« argument, among other issues, is the increased static stability induced by mid-troposphere 

* warming and surface cooling. Similar findings are also obtained by Lau and Kim [2007] 

* on a seasonal time-scale. Reale et al. [2009b], hereafter RA09b, reiving upon a global 
data assimilation and modeling effort, and the production of high quality analyses, also 

58 suggest that high dust content is associated with a. thermal dipole: relatively warm at 
» 600-700 hPa and cooler at about 900 hPa or below. RAOOb findings are supported by 
w the fact that the finite- volume dynamics of the NASA GEOS-5 is particular suitable to 
ii maintain fine thermal features avoiding unrealistic dispersions. 
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Recent climate sensitivity experiments using general circnlaf ion models have shown, t hat 
radiative effects of Saharan dust tend to draw the Atlantic ITCZ northward towards the 
southern edge of the dust layer, and induce large-scale north-south and. east-west oriented 
divergent circulations. that affect the African Easterly Waives and rainfall variability of the 
West African monsoon and climate of the Atlantic l.nu i t. ah, 2009: Kim ef a].: 201 Oj. 
Wilcox et al. [2010] have shown from observational perspective that the temperature 
structure associated with Saharan air outbreaks is consistent with previous work, and that 
one of the effects of dust outbreaks is a northward shift, of the filter Tropical Convergence 
Zone (ITCZ). Also within the operational forecasting modeling framework, evidence has 
been provided by Tompkins et al. [2005] on that an improved representation of aerosol (i.e. 
seasonally varying instead of fixed) leads to a. more correct, representation of the African 
Easterly Jet (namely, a northward and upward shift, in agreement with observations). 

In this work, the use of interactive aerosol in a global modeling framework provides 
results fully consistent with Tan et al. [2009], Kim et al. [2010]', Wilcox et al. [2010] 
and confirms that the direction started by the pioneering work by Tompkins et al. [2005] 
needs to be further pursued. 


2. The Model and Data Assimilation System 

The global data assimilation and forecasting system used is the NASA GEOS-5. which 
combines the Griclpoint Statistical Interpolation (GSI) analysis algorithm developed by 
the National Centers for Environmental Predictions (NCEP) (e.g. Wu et ah, [2002]). and 
modified by the NASA Global Modeling Assimilation Office (GMAO). with the NASA 
atmospheric global forecast model, is documented in Rieuecker et al.. [2008]. The fore- 
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cast model slum*:- the same dynamical core '[.in. 200 -If with the so-called finite- volume 


General ( Irculntion Model (fvGCM), also known, as. tin*. GEOS- 1. Inn contains a different 


set of physical pannnetrixal ions, partly developed by the NASA GM AO. The Aerosol 
interactive component is a new feature which is currently implemented also in the opera- 
tional version of the CEOS-5 and is a further step with respect, to the previous transport 
model described in Colareo et al. [201 0J. which merged the Goddard Chemistry, Aerosol, 
Radiation and Transport Model (GOCART) with the. GEOS-4. The version described 
by Colareo et al. [2010] allows the treatment of dust, sea. salt, carbonaceous and sulfate 
aerosols, producing realistic aerosol distributions (validated against ground- and satellite- 
based measurements) which are consistent with the dynamics and the meteorological 
fields produced by the GEOS-4. In the version used in this study, aerosol transport, pro- 
cesses (adveetion, diffusion, convection) are provided by the host GEOS-5 model, with the 
aerosol specific processes (emission, deposition, simplified sulfate chemistry) as in Colareo 
et al. [2010]. The direct aerosol radiative effect, is included in the GEOS-5, while indirect 
effects on cloud/precipitation processes are not included at present. 


3. The Experiments 

Four sets of 30 5-day forecast, experiments at two horizontal resolutions (0.5° x 0.67" 
and 0.25" x 0.33", both with 72 vertical levels), are performed. They are initialized 
daily at OOz, starting at OOz 15 August 2006. from the same set of analyses used in 
RA09b. These were produced by assimilating all. conventional and satellite observat ions 
used operationally at that rime, but with a. much denser coverage from the Atmospheric 
Infrared Sounder (AIRS), obtained by ingesting cloudy retrievals instead of dear-sky 
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radiances. As discussed in Reale ef al. 1 2009a] and in HA091), AIRS-dorived information 
iu partly cloudy regions improves the analyses jn the tropics. However, the experiments 
described in this article are characterized by the inclusion or exclusion of interactive 
aerosol, and are named N()A rx (no aerosol) and f A A tx (interactive aerosol) with the 
suffix indicating the 0.50 l! or 0.25" forecasts respectively. 


4. Results 

As a preliminary metric to assess the impact of interactive aerosol, the anomaly corre- 
lation for all sets of 30 5-day forecasts is computed against operational NCEP analyses. 
The impact is virtually zero: the insertion of interactive aerosol does not affect the globally 
averaged skill, of the system. In PS I the anomaly correlations of .WA R,, and 1 AA r> o are 
compared; the 0.25° forecasts provide identical results ( not shown). The lack of global 
impact of the interactive aerosol is understandable because the production of dust is con- 
fined to few areas of the world, and does not occur al all times even in the dust-producing 
regions. However, it is important to notice that the forecast skill does not degrade as a 
consequence of the insertion of aerosols, which represents a major change in the model's 
physics. 

At the same time, it is crucial to emphasize that, for states of the atmospheres associated 
with specific weather systems, at times in which dust production is strong, and selected 
regions are analyzed, the impact- of aerosol insertion is very large and dynamically relevant. 
In this work we focus on the impact on the African Easterly Jet (AEJ) which flows, 
with a predominantly easterly component, at about 6-50 hPa and about- 15" — 20° A . and 
which is partly controlled bv horizontal temperature gradients over the Sahelian region. 
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To detect and understand the impact of interactive aerosols is not an easy or obvious 
task, because of the intrinsically noisy nature of the radiative forcing associated with 
dust, t he overlapping with strong thermal signal due to the diurnal cycle, and finally the 
inhomogeneous distribution of dust, rapidly varying in space and time. In this work, wo 
aim to .show an example of a dear signal associated to a major dust outbreak:. We select 
the strong outbreak observed, during the Special Observing Phase (SOP-3) of NAMMA, 
moving from Africa- to the Atlantic between 25 and 28 August. RAOOb suggest that 
this (went produced a temperature dipole (cool at about, 900 hPa and warm at about 600 
hPa), whose effect propagated well into the At lantic. However, R.A09b did not have yet; the 
capability of interactive aerosol and is based on a comparison between the MODIS AODs, 
and on. the temperature structure in the GEOS-, 5 analyses and integrations [Reale and 
Lau 2010]. In this work we rediscuss the presence of a t hermal dipole possibly associated 
with dust, in lieu of this new NASA tool. 

As a first step, it is important to verify that the above-referred dust, outbreak is rep- 
resented well in the simulations. Figure 1 com part's the outbreak as represented through 
the total mass column dust obtained from the 36-hour GEOS-5 forecast initialized at OOz 
24 August, and a, composition of AOTs obtained from MODIS instruments on board both 
Aqua and Terra satellites, corresponding to the verification time of 12z 2-5 August. The 
Figure shows that the GEOS-5 reproduces a strong dust outbreak intercepting the African 
coastline at about 20° — 30" A 7 , recurving southward at about 200!’ and northwestward 
at. about. 23° IE — 35°1 1 
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Figure FS2 -shows the same outbreak, as raptured by a 30-hour forecast initialized on 
the following day, at 00? 25 August 2006. By comparing the observed AOD from 12?. 26 
>.v August with the forecast, it can be seen that GEOS-5 again captures well the westward 
ix. propagat ion of the dust edge, and even produces a very thin low-dust: channel between the 
two major high dust regions centered at about. 15" A . 30 — 40" U’ and 25" A'. 20 — 30' ’ll . 
m respectively. On these grounds, it can be safely stated that there is a major dust out break 
.I? iti. t!.io GEOS-5, and that its timing and scale correspond well with the MODIS data, 
in The second step is to investigate the impact of this dust outbreak on the temperature 
im structure, as represented by the GEOS-5. Figure 2 shows the temperature impact of 
us aerosols at longitude 10'TF, from the same simulations in Fig. 1. An evident strong 
im thermal dipole, very similar to the one discussed in RA09b and previously by Dunion and 
ui Widen [2004] can be observed: a healing of about V’C - \.5°C between 850 and 500 liPa, 
i38 and a cooling in the lower, near-surface, levels. A very similar pattern in the IAA& minus 
me NOAM temperature structure across the same latitude range and at longitude KBIT can 
no be seen also in FS 3, where the corresponding aerosol impact on temperature is computed 
hi from the 36 hour forecast initialized on the following day (00z 25 August } and verified at 
i« 12? 26 August. The simulations at 0.50" provide consistent results (not shown) but with 
m diluted values and smoother gradients. As discussed in RA09b and Reale and Lau (2010). 
us the resolution of 0.25" is better than 0.50" to investigate the fine structure of the AEJ. 
i« For this reason, we show only the IAA& in the following pictures. 
t« Because the AEJ and the wind over the region are strongly controlled by lower-level 
i« temperature gradients, a change in the thermal structure must have an impact on the 
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iW AEJ ii ih 1 more generally on ( he ciix ul.at.ioa. Figure 3 shows the vertical profiles extracted 
n<> from a 108 hour forecast in t he GEOS-5 integrations i A ( > 1 am! / 1.1 :■ and compared 
with the profiles of temperature and zonal, wind component, obtained from the vertical 
iii ground-based soundings from the Cape Verde Islands, at i 23.,Y'il ' (radiosonde 

i« data available online at: http: / /namn ia. i n.sh ■. n asa. gov ) , Tlie most: evident impact on 
■* temperature is a warming up to l°f ’ between 900 and 700 hPa, and a slight cooling from 
«♦ the near-surface levels up to 900 hPa. The zonal wind panel shows that the verifying 
iw analysis and the observational profile match very well, but the 108 hour NOA-^r, forecast 
is. has a negative bias with respect to the observations, of about 4ms " 1 . from the surface to 
t« SOOhPa. 1 he impact ol the aerosol is very evklent in the lower levels, reducing the bias of 
ir,9 about 40%. When no dust production is present, the impact on temperature and zonal 
i» wind is negligible (not shown}. Also, above OODhPa. the impact, of the aerosol is virtually 
iso zero (not shown). 

im The radiative effect of dust is controlled by diurnal cycle and the distribution of aerosol 
>« is very noisy; therefore, to obtain a clearer signal on the wind field and of the AEJ 
m structure, a 5-day average across the forecast initialized on August 25th is performed. 
i« There is an evident northward and upward shift in the AEJ structure, in the zonal wind 
and temperature fields at 10' if (FS 4) consequent to this thermal anomaly, in agreement 
nk. with Tompkins et al. (2005!, and expected from the thermal wind relationship. We 
I*, ■ see a similar thermal anomaly at 1GPE and a corresponding northward shift in the AEJ 
«« (FS 5). The shift is also very evident on any longitude between 20*11’ and 20*1? (not 
i6» shown) and, to a. lesser extent, also in the lAA r> a experiments (not shown). The scarcity 
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of observations from the arm (which prevents from providing a dew depiction of the 
AE.J from observations only) does not allow the validation of this northward shift in a 
more rigorous war. To compare the AEJ with reanalyzes is also not meaningful since 
no re; in a lysis is produced from a .system containing realistic, interactive aerosol at this 
time. However, these experiments demonstrate that the effect: of interact ive aerosol in the 
model, without degrading the global forecasting skill, alter .significantly the circulation 
over the dust-affected region and provide an improved vertical profile at one validation 
point, 


5. Discussion 

The introduction of interactive aerosol, varying in agreement with the meteorology, is a 
serious challenge requiring consideration, in this work it is shown that the usual metrics 
to asses improvement in forecasting skill, such as 500 hPa anomaly correlation plots, are 
not. suitable to evaluate the impact of the insertion of interactive- aerosol treatment. On 
the contrary, for areas and events affected by dust aerosols, if is important to design mean- 
ingful, event-focused, metrics. In this work a major dust outbreak observed during the 
NAMMA campaign in August 2006 is investigated with the aid of the GEOS-5 forecasting 
system inclusive of interactive aerosol component. It is first shown that the forecast, of 
dust total column mass qualitatively agrees with the MODIS merged AOTs from both 
Aqua and Terra, and then that the effect of the aerosol on the temperature structure is 
completely consistent with previous observational and modeling works, producing a ther- 
mal dipole (warmer at the aerosol level, cooler below}, hi particular, if is shown that the 
effect: of this thermal anomaly projects into the dynamics and produces an improved 108 
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hour foi r-f.!-! of temperature and wind profiles, validated against mi existing sounding. 
A 5-day average across the forecast, to dampen the noisy structure of the dust, -induces 
signal, shows that the net effect, of the aerosol simulation is a northward (and. to a lesser 
extent , upward) shift of the AEJ. in agreement with other modeling studies. FS i and 5 
also suggest a nort hward displacement of the westerly low-level mousoouaJ flows in agree- 
ment with Wilcox et al. [2010] and related modeling results [e.g. Lmi et al 2009]. In 
this work we have also performed the same analysis for another dust outbreak occurred 
around 10-13 September: wind and temperature anomalies of the same sign are obtained, 
albeit of a weaker magnitude, being the outbreak also substantially weaker (not shown). 


6. Concluding Remarks 

In this article we show that the use of interactive aerosol in t he GEOS-5 data assimila- 
tion and forecasting system does not, impact the global skill but; substantially changes the 
thermal structure over northwestern tropical Africa, during an observed dust outbreak, 
affecting t he circulation and particularly the AEJ representation. We validate the changes 
in temperature and wind with a. vertical sounding, showing f hat the insertion of interactive 
aerosol improves the forecast up to day five for this particular event. However, the limi- 
tation of this work is to rely upon modeled aerosol production and distribution. The next 
step, namely real-time use of satellite-derived aerosol information in global, models, is a 
very rec ent development of the most advanced data assimilation systems. In this regard, 
the NASA CM AO is currently attempting to assimilate real-time AOTs from MODIS. 
which may represent the future direction for this kind of research. 
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Figure 1 . 36-hour forecast (initialized at OOz 24 Aug 2006} of total mass column aerosol from 


the GEOS-5, and MOD IS optical depth for 12z 25 Aug 2006. 
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Figure 2. Vertical Meridional cross-section of temperature ( oC ) at 10° W in the GEOS-5 
A* O/l 25 simulation for 12z 25 Aug 2006, (36 hour forecast, initialized at OOz 24 Aug). Shaded: 
/ A A 25 minus NOA^ difference. 
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Figure 3. Vertical profiles of Temperature (°C. left panel) and zonal wind u (ms" 1 , right 
panel) at Cape Verde 14.92°A r , 23.49‘TS' in the GEOS-5 NOA-a, (green) simulation for 12z 29 
Aug 2006, (108 hour forecast, initialized at OOz 25 Aug), verifying analyses (black), compared 
with observed soundings (purple) at the Cape Verde, from SOP-3 campaign. In the left panel the 
orange line represents the MA25 minus NOA% 5 difference; in the right panel, the actual IAA& 
u, component of the wind. 
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